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VLBI observations of relativistic outflows (jets) in galactic nuclei, complemented with detailed 
studies made in other spectral domains, have become an effective tool for investigating the physics 
of nuclear regions in galaxies. High-resolution radio observations access directly the regions 
where the jets are formed, and trace their evolution and interaction with the nuclear environment. 
The emission properties, dynamics, and evolution of jets in AGN are intimately connected to the 
characteristics of the central supermassive black hole, accretion disk, and broad-line region (BLR) 
in active galaxies. Large VLBI surveys (15 GHz VLBA survey, MOJAVE) and dedicated mon- 
itoring programmes follow systematically the evolution of several hundreds of relativistic jets. 
These observations, combined with optical and X-ray studies, yield arguably the most detailed 
picture of the galactic nuclei. Recent results from studies of the nuclear regions in several active 
galaxies with prominent outflows are reviewed in this contribution. 
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Jets and nuclear regions in active galaxies 



1. Introduction 



Substantial progress achieved during the past decade in studies of active galactic nuclei (see [51] 
for a review of recent results) has brought an increasingly wider recognition of the ubiquity of 
relativistic outflows (jets) in galactic nuclei [|l7|, |84| |, turning them into an effective probe of nu- 



clear regions in galaxies [46]. Emission properties, dynamics, and evolution of an extragalactic 
jet are intimately connected to the characteristics of the supermassive black hole, accretion disk, 



and broad-line region (BLR) in the nucleus of the host galaxy [51]. The jet continuum emission 
is dominated by non-thermal synchrotron and inverse-Compton radiation [f78|]. The synchrotron 
mechanism plays a more prominent role in the radio domain and the properties of the emitting 



material can be assessed using the turnover point in the synchrotron spectrum Q45|], synchrotron 
self-absorption [44], and free-free absorption in the plasma [31, [Si] ]. 



High-resolution radio observations access directly the regions where the jets are formed [3C] 



and trace their evolution and interaction with the nuclear environment [62]. Evolution of compact 
radio emission from several hundreds of extragalactic relativistic jets is now systematically studied 
with dedicated monitoring programmes and large surveys using very long baseline interferometry 



(VLBI) such as the 15 GHz VLB A 1 survey Q and MOJAVE [Q. These studies, combined 
with optical and X-ray studies, yield arguably the most detailed picture of the galactic nuclei [pq]. 
Presented below is a brief summary of recent results in this field, outlining the relation between jets, 
supermassive black holes, accretion disks, and BLR in prominent active galactic nuclei (AGN). In 
this respect, this review is complementary to other recent reviews [[H], |8|] focused on formation 
and propagation of extragalactic relativistic jets. 



2. Anatomy of jets 



Jets in active galaxies are formed in the immediate vicinity of the central black hole [ |l0| ] and 
they interact with every major constituent of AGN (see Figure [l] and Table [l]). The jets carry away 
a fraction of the angular momentum and energy stored in the accretion flow p3] or corona (in low 



luminosity AGN Q61TJ) and in the rotation of the central black hole [ |37| , p8[ [72f J . 

The production of highly-relativistic outflows requires a large fraction of the energy to be 



converted to Poynting flux in the very central region [73]. The efficiency of this process may 



depend on the spin of the central black hole [59]. The collimation of such a jet requires either 



a large-scale poloidal magnetic field threading the disk [ |75[ ] or a slower and more massive MHD 
outflow launched at larger disk radii by centrifugal forces [9[ 77]. A two-zone version of such a 



hybrid outflow model is known as the "two-fluid model" [|74[]. 

At distances of ~ 10 3 /? g (R g = GM/c 2 is the gravitational radius of a black hole), the jets 
become visible in the radio regime, which makes high-resolution VLBI observations a tool of 



choice for probing directly the physical conditions in AGN on such small scales [30, 39]. Recent 
studies indicate that at distances of 10 3 -10 5 /? g (< lpc) the jets are likely to be dominated by 



pure electromagnetic processes such as Poynting flux [73] or have both MHD and electrodynamic 
components [^OQ. The flowing plasma is likely to be dominated by electron-positron pairs [82, 
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Figure 1: Basic sketch of a radio-loud AGN featuring a relativistic jet and its environment. The scale 
is logarithmic beyond 10/? g and vertical dimensions are not drawn to scale. Regions of the jet mainly 
contributing to various parts of the broadband spectrum are indicated. Adapted from 



Table 1: Characteristic scales in the nuclear regions in active galaxies 
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h 


^Gpc 




^orb 




[*g] 


[pc] 


[mas] 


[yr] 


[yr] 


Event horizon: 


1-2 


i(r 5 


5 x 10~ 6 


0.0001 


0.001 


Ergosphere: 


1-2 


lo -5 


5 x 10~ 6 


0.0001 


0.001 


Corona: 


lO'-lO 2 


10- 4 -10- 3 


5 x lo" 4 


0.001-0.01 


0.2-0.5 


Accretion disk: 


lO'-lO 3 


10- 4 -10- 2 


0.005 


0.001-0.1 


0.2-15 


Jet formation: 


>10 2 


>io- 3 


>5 x lO- 4 


>0.01 


>0.5 


Jet visible in the radio: 


>10 3 


>io- 2 


>0.005 


>0.1 


>15 


Broad-line region: 


10 2 -10 5 


10 3 -l 


0.05 


0.01-10 


0.5-15000 


Molecular torus: 


>10 5 


>1 


>0.5 


>10 


> 15000 


Narrow-line region: 


>10 6 


>10 


>5 


>100 


>500000 



Column designation: I - dimensionless scale in units of the gravitational radius, GM/c 2 ; 1% - corresponding 
linear scale, for a black hole with a mass of 5 x 1O 8 M ; 0Gpc - corresponding largest angular scale at 1 Gpc 
distance; x c - rest frame light crossing time; T or b - rest frame orbital period, for a circular Keplerian orbit. 
Adapted from [|1} 



23], although a dynamically significant proton component cannot be completely ruled out at the 



moment [14]. 

The magnetic field is believed to play an important role in accelerating and collimating ex- 
tragalactic jets on parsec scales [80]. How far the magnetic field dominated region extends in 
extragalactic jet is still a matter of debate [[73TJ. Nevertheless, it is possible to identify three dis- 
tinct regions with different physical mechanisms dominating the observed properties of the jets: 
1) ultracompact jets where collimation and acceleration of the flow occurs, 2) par sec-scale flows 
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dominated by relativistic shocks, and 3) large-scale jets where plasma instabilities dominate the 
flow. 




Figure 2: Left panel: Relative changes of the Doppler factor and magnetic field in the VLBI core in 3C 345, 
obtained by applying Konigl jet model to the measured flux density and frequency of the synchrotron peak 
in the spectrum. All quantities are normalized to their respective values at the first epoch, to — 1981 .5. Open 
circles denote the particle density required for maintaining a constant Doppler factor. The dotted line shows 
the particle density, as represented by five exponential flares. The resulting Doppler factor (open squares) and 
magnetic field (filled triangles) are also shown, with lines representing linear fits to the respective quantity. 
The light curve of the total flux density at 22 GHz is plotted for comparison, scaled down by a factor of 100. 
Reproduced from [p^. Right panel: Values of apparent transverse speed, j3 app , and apparent luminosity, L 
plotted for the fastest component of 129 radio sources in the 15 GHz VLB A survey [[33]]. Colours indicate 
different host galaxy types. The aspect curves are the loci of (j3 app ,£) for sources with Lorentz factors Tj = 30 
and 20, andL = 1 x lO^WHz -1 . Reproduced from [@. 



2.1 Ultracompact jets 



Ultracompact jets observed down to sub-parsec scales typically show strongly variable but 
weakly polarized emission (possibly due to limited resolution of the observations). Compelling 
evidence exists for acceleration [§] and collimation [60, |39|, B0] of the flows on these scales, which 



is most likely driven by the magnetic field []8Q|]. The ultracompact outflows are probably dominated 



by electromagnetic processes [|60[ [73J] and they become visible in the radio regime (identified as 
compact "cores" of jets ) at the point where the jet becomes optically thin for radio emission [44, 
^9h. The ultracompact jets do not appear to have strong shocks [^] and their basic properties 
are successfully described by quasi-stationary flows [[35]]. The evolution and variability of ultra- 
compact jets (Figure Q left panel) can be explained by smooth changes in particle density of the 
flowing plasma, associated with the nuclear flares in the central engine J49|]. Intrinsic brightness 



temperatures of the ultracompact jets are estimated to be (1-5) x 10 11 K [54], implying that the 



energy losses are dominated by the inverse-Compton process [|32|]. 

Quasi-periodic variability of the radio emission from the ultracompact jets is most likely re- 
lated to instabilities and non-stationary processes in the accretion disks around central black holes 
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Figure 3: Left panel: Distribution of the synchrotron turnover frequency, V m , in the jet in 3C273. The 
distribution is obtained from multifrequency VLBI observations and shows three regions of increased v m : 
one is in the "core" of the jet and the other two are coincident with two bright superluminal features in the 
jet (denoted Jl and J2). Right panel: Profiles of the turnover frequency (blue) and magnetic field strength 
(red) taken in the jet, along the dashed line shown in the left panel. The low magnetic field in the core 
region indicates that the field is tangled or intrinsically weak there. The magnetic field profile shows strong 
spikes at the locations of the features Jl and J2, implying that these are most likely due to magnetic field 
compression in strong shocks propagating in the jet. Adapted from 



in AGN [|26|, fig]. Alternative explanations involve binary black hole systems in which flares are 
caused by passages of the secondary through the accretion disk around the primary [27, ^2]. These 



models, however, require very tight binary systems, with orbits of the secondary lying well within 
10 3 a Schwarzschild radii of the primary (between 20 and 100 Schwarzschild radii, in the celebrated 



case of OJ287 Q42|]). This poses inevitable problems for maintaining an accretion disk around the 
primary (for massive secondaries [^]) or rapid alignment of the secondary with the plane of the 
accretion disk (for small secondaries [^]). Discrepancy between the predicted and actual epoch of 
the latest outburst in OJ 287 [[7^] indicates further that the observed behaviour is not easy to be re- 
produced by a binary black hole scenario and it is indeed more likely to result from a quasi-periodic 



process in the disk, similarly to the flaring activity observed in 3C 345 [ ]48[ ] . 

2.2 Parsec-scale flows: shocks and instabilities 

Parsec-scale outflows are characterized by pronounced curvature of trajectories of superlumi- 
nal components [^53|, &9|], rapid changes of velocity and flux density, and predominantly transverse 
magnetic field ||29|]. Statistical studies of speed and brightness temperature distributions observed 
in the superluminal features propagating on parsec scales indicate that the jet population has an en- 
velope Lorentz factor of 30 and an unbeamed luminosity of ~ 1 x 10 25 W Hz -1 [ flr^ ] (right panel 
of Figure with the brightness temperature of the emitting plasma reaching ~ 5 x 10 10 K [^4|], 
which is close to the equipartition limit [|70|]. 

Relativistic shocks are expected to be prominent on these scales, which is manifested by strong 
polarization | f7T| ] and rapid evolution of the turnover frequency of synchrotron emission [52]. Map- 
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Figure 4: Left panel: The observed trajectory of C5 in the plane of the sky. The dotted line represents the 
combined polynomial fits to the component's x and y offsets from the core. Open circles mark the locations 
on the trajectory which are equally spaced in time at an interval of 1 year. Right panel: Evolution of the 
synchrotron turnover changes in C5. The dotted line shows an evolutionary track consistent with different 
stages of the evolution of the shock (Compton-loss, synchrotron-loss, and adiabatic-loss stage). The dot- 
dashed line shows how the original fit for C5 must be changed to satisfy also the observed trajectory of 
the component (shown in the left panel). The shock description can only be applied at early stages of the 
component evolution, indicating that the shock should dissipate rapidly, on scales of < 100 pc. Adapted 
from [|9|. 



ping the turnover frequency distribution provides also a sensitive diagnostic of shocks and plasma 
instabilities in relativistic flows [45]. Shocks are particularly evident in the profiles of magnetic 
field obtained from the turnover frequency images (see Figure [3]). 

Complex evolution of shocked regions is revealed in observations [|I^, |29|, ^9|] and numerical 
simulations [J3j] of parsec-scale outflows. However, the shocks are shown to dissipate rapidly [ |49| ] 
and shock-dominated regions are not likely to extend beyond ~ 100 pc. This can be exemplified by 
the kinematic and spectral evolution observed in the jet component C5 in 3C 345 (Figure |j). The 
shock model can reproduce the observed spectral and kinematic changes only during the initial 
stages of the component evolution. This implies that strong shocks should dissipate rapidly in the 
jets, on scales of < 100 pc [49]. 

Observed morphologies of parsec-scale jets and trajectories of superluminal features propa- 
gating in the jets are often described in terms of a helical geometry j76|], with helicity arising from 



some periodic process in the nucleus. Jet precession, both in single [55] and binary [13] black hole 
systems, have been commonly sought to be responsible for the observed helicity on parsec scales, 
but it may lead to rather non-physical parameters of the nuclear region (see [48] for discussion). 
It is more likely that the jet precession in AGN acts on time scales of > 10 4 years and manifests 
itself in kiloparsec-scale structures [ 21 ] where it also becomes affected by the motion of the host 
galaxy [83]. Determining the precession parameters from observations of parsec-scale structures 
requires detailed description of relativistic effects and source geometry. This yields much longer 
orbital and precession periods [48] than those obtained from associating quasi-periodic changes ob- 
served on these scales directly with precession [|13|]. On parsec scales, the observed helical patterns 
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Figure 5: Top: HST image of the jet in M 87. Bottom: Model based on linear description of Kelvin- 
Helmholtz instability developing in the jet. Line-of-sight synchrotron intensity image and contours at a jet 
viewing angle of 40° including all light travel time effects and time delays. The beginning point is about 2" 
out from the nucleus with the end point at about 22". The spatial twist is such that the twisted flow is more 
towards the observer at the bottom of the jet and is more away at the top of the jet. Thus, knots A and C 
would be Doppler de-boosted and knot B would be Doppler boosted. Reproduced from [p3||. 



may be produced not exclusively by precession but also by other processes, including Kelvin- 
Helmholtz instability [ p2| ] and rotation of the flow [fTl]]. The flow rotation is particularly relevant 
for explaining periodic changes of the "ejection angle" (flow direction in the immediate vicinity of 
the VLBI core on sub-milliarcsecond scales) reported in several prominent objects [jlj, ^3[]. 

2.3 Large-scale jets 

On scales larger than ~ 100 pc, instabilities (most importantly, Kelvin-Helmholtz instability) 
determine at large the observed structure and dynamics of extragalactic jets [50, |3|, 59]. The ellip- 
tical mode of the instability is responsible for appearance of thread-like features in the jet interior, 
while overall oscillations of the jet ridge line are caused by the helical surface mode. Success- 
ful attempts have been made to represent the observed brightness distribution of radio emission 
on these scales, using linear perturbation theory of Kelvin-Helmholtz instability [ p3| ] and a spine- 
sheath (analogous to the two-fluid) description of relativistic flows [|TJ, ffj]|. Figure |5| illustrates 
application of linear Kelvin-Helmholtz model to reproducing the radio brightness distribution in 
the jet of M87 Non-linear evolution of Kelvin-Helmholtz instability [ ]66[ 57] and stratifica- 
tion of the flow [ ]68| ] play important roles in large-scale jets. Similarly to stellar jets, rotation of the 



flow is expected to be important for extragalactic jets Q18[], but observational evidence remains very 
limited on this issue. 



3. Jets and nuclear regions in AGN 

A number of recent studies have used jets to probe physical conditions in the central regions 
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Figure 6: Left panel: Magnetic field distribution in the jet in 3C 345. Squares show the magnetic field in 
the compact jet derived from the frequency-dependent shift of the core. Circles are the homogeneous syn- 
chrotron model estimates of magnetic field in the extended jet components. Triangles show the characteristic 
magnetic field values from a model of magnetized accretion disk. Diamonds are the theoretical estimates 
for the dipole magnetic field around a supermassive rotating black hole. Right panel: Opacity in the jet 
in 3C 309.1. Circles are the measured values of the core shift power index k T in 3C 309.1 at different fre- 
quencies; solid line shows changes of k r due to pressure gradients in the BLR clouds supported by thermal 
pressure and maintaining a mass distribution with spherically symmetrical gravitational potential. The cloud 
region extends up to 400 r s (r s refers to the distance at which the jet becomes supersonic). The equipartition 
regime is approached at the outer boundary of the cloud region, with k r = 1 . Significant deviations from 
the equipartition are seen on smaller scales, resulting in stronger self-absorption in the inner parts of the jet. 
Reproduced from |Q] . 



of AGN. These studies have probed physical conditions in the compact relativistic flows, proper- 
ties of atomic and molecular material in circumnuclear regions of AGN, and connection between 
relativistic outflows, accretion disks, and BLRs. 



3.1 Jets and nuclear regions of AGN 

Synchrotron self-absorption and external absorption in the ultracompact jets (VLBI "cores") 
can be used effectively for determining the properties of the flow itself and its environment [H]. 
Absolute position of the core, r c , varies with the observing frequency, v, so that r c oc y _1 /*% [p5|]. 
If the core is self-absorbed and in equipartition, the power index k r = 1 Changes of the core 
position measured between three or more frequencies can be used for determining the value of 
k x and estimating the strength of the magnetic field, B cor e> in the nuclear region and the off- 
set, /?core> of the observed core positions from the true base of the jet (see Figure ^; left panel). 
The combination of these gives an estimate for the mass of the central black hole Mbh ~ 7 x 

10 9 M (B C ore/G) l ' 2 (/^ore/PC) 3 / 2 . 

Core shift measurements provide estimates of the total (kinetic + magnetic field) power, the 
synchrotron luminosity, and the maximum brightness temperature, 7b, max in the jets. The ratio 
of particle energy to magnetic field energy can also be estimated, from the derived 7b,max- This 
enables testing the original Konigl model [35] and several of its later modifications (e.g., [25, 8p). 
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The known distance from the nucleus to the jet origin will also enable constraining the self-similar 
jet model [ |56| ] and the particle-cascade model [Q]. 

Recent studies of free-free absorption in AGN indicate the presence of dense, ionized cir- 
cumnuclear material with T e « 10 4 K distributed within a fraction of a parsec of the central nu- 



cleus [|44j, |81p . Properties of the circumnuclear material can also be studied using the variability 
of the power index k T with frequency. This variability results from pressure and density gradients 
or absorption in the surrounding medium, most likely associated with BLR. Changes of k T with 
frequency, if measured with required precision, can be used to estimate the size, particle density, 
and temperature of the absorbing material surrounding the jets (see Figure ||; right panel). Esti- 
mates of the black hole mass and the BLR size obtained from the core-shift measurements can be 
compared with the respective estimates obtained from the reverberation mapping and applications 
of the Mbh-cr* relation. 

3.2 Atomic and molecular absorption 

Opacity and absoiption in the nuclear regions of AGN have been probed effectively using the 



non-thermal continuum emission as a background source Q46[]. Absorption due to several atomic 
and molecular species (most notably due to H i, CO, OH, and HCO + ) has been detected in a number 
of extragalactic objects. OH absorption has been used to probe the conditions in warm neutral 



gas [ 201 , 34 ], and CO and Hi absorption have been used to study the molecular tori []15J, |65|] at a 
linear resolution often smaller than a parsec [|2|]. These observations have revealed the presence of 
neutral gas in a molecular torus in NGC4151 and in a rotating outflow surrounding the relativistic 



jet in 1946+708 [64] 



3.3 Jet-disk and jet-BLR connections 



Connection between accretion disks and relativistic outflows [24] has been explored using 



correlations between variability of X-ray emission produced in the inner regions of accretion disks 



and ejections of relativistic plasma into the flow [57]. The jets can also play a role in the generation 



of broad emission lines in AGN. The beamed continuum emission from relativistic jet plasma can 
illuminate atomic material moving in a sub-relativistic outflow from the nucleus, producing broad 
line emission in a conically shaped region located at a significant distance above the accretion 
disk Magnetically confined outflows can also contain information about the dynamic evolution 
of the central engine, for instance that of a binary black hole system p8|]. This approach can be 
used for explaining, within a single framework, the observed optical variability, kinematics, and 
flux density changes of superluminal features embedded in radio jets. 



4. Conclusion 

Extragalactic jets are an excellent laboratory for studying physics of relativistic outflows and 
probing conditions in the central regions of active galaxies. Recent studies of extragalactic jets 
show that they are formed in the immediate vicinity of central black holes in galaxies and carry 
away a substantial fraction of the angular momentum and energy stored in the accretion flow and 
rotation of the black hole. The jets are most likely collimated and accelerated by electromag- 
netic fields. Relativistic shocks are present in the flows on small scales, but dissipate on scales 
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of < 100 pc. Plasma instabilities dominate the flows on larger scales. Convincing observational 
evidence exists, connecting ejections of material into the flow with instabilities in the inner accre- 
tion disks. In radio-loud objects, continuum emission from the jets may also drive broad emission 
lines generated in sub-relativistic outflows surrounding the jets. Magnetically confined outflows 
may preserve information about the dynamics state of the central region, allowing detailed inves- 
tigations of jet precession and binary black hole evolution to be made. This makes studies of 
extragalactic jets a powerful tool for addressing the general questions of physics and evolution of 
nuclear activity in galaxies. 
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